Superconductivity in non-centrosymmetric BiPd system 
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In this work, we establish the bulk superconductivity of a high quality sample of monoclinic BiPd (a-BiPd, 
space group P2[) below 3.87 K by studying its electrical resistivity, magnetic susceptibility and heat capacity. 
We show that it is clean type-II superconductor with moderate electron-phonon coupling and determine its 
superconducitng and normal state parameters. Although o--BiPd is a noncentrosymmetric superconductor with 
large electronic heat capacity (therefore, large y), the effect of spin-orbit splitting of the electronic bands at the 
Fermi level is small. This makes little influence on the superconducting properties of a-BiPd. 
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I. INTRODUCTION 



Ever since the discovery of the non-centrosymmetric heavy 
fermion superconductor CePt^Si 1 , there is widespread re- 
search activity to understand the nature of superconductiv- 
ity in such unconventional superconductors. The term "non- 
centrosymmetric" characterizes the symmetry of a crystal lat- 
tice without inversion center. In such materials the standard 
classification in even-parity spin-singlet and odd-parity spin 
triplet superconducting phases is obsolete, because the elec- 
trons are exposed to antisymmetric spin-orbit coupling, e.g. 
Rashba-type of coupling 2 which arises due to electric field 
gradient in the crystal which has no inversion symmetry. An 
inherent feature is then the mixing of spin-singlet and spin- 
triplet Cooper pairing channels which are otherwise distin- 
guished by parity. This mixing of pairing states is expected 
to cause a two-component order parameter. New forms of 
pairing appear giving rise to unusual temperature/field depen- 
dence of the superconducting parameters. During recent years 
we have witnessed rapid developments on the side of experi- 
ments and synthesis of novel materials as well as in the theo- 
retical understanding of this type of superconductors. Indeed 
many new materials have been found, in particular, among 
the heavy fermions for which unconventional Cooper pair- 
ing is expected. Exotic mixed-state phases (vortex matter) 
have been predicted in theory and are under experimental in- 
vestigations. For example in the case of CePfjSi where su- 
perconductivity occurs at ambient pressure, whereas, transi- 
tion to such a state happens in Ulr— , CeRhSiji and CelrSi^ 
only under pressure. However, the study of superconductiv- 
ity in non-centrosymmetric materials which do not exhibit 
heavy fermion features is also important since it avoids ad- 
ditional complication that arises due strong f-electron corre- 
lations. Discovery of such materials also continues to in- 
crease, starting from binary carbides (R 2 C3_,5 with R=La or 
Y)£, Cd 2 Re 2 Ori , Li 2 (Pd,Pt) 3 B^, Mg 2 Al^ and the recently 
found BaPtSijii. However, many of them (except Li 2 Pt3B) 
exhibit conventional BCS-like superconductivity due to small 
spin-orbit scattering. One of the common features in these 
compounds which show conventional superconductivity with 



small spin-orbit scattering is the absence of the large density 
of states at the Fermi level. Hence, it will be of interest to 
study a superconducting material which have conduction elec- 
trons with high density of states at the Fermi level (not from 
high f-electron correlations) but has no inversion symmetry. 
In this work we report our investigations in one such mate- 
rial, namely c-BiPd, that shows bulk superconductivity be- 
low 3.87 K via, resistivity, magnetization and heat capacity 
studies. ff-BiPd has no inversion symmetry in its monoclinic 
crystal structure (space group P2i). One of the important fea- 
tures of the structure is that the presence of short Pd-Pd bonds 
(shorter than those present in pure Pd metal). This could re- 
sults in conduction electrons with large density of states at the 
Fermi level which is suggested in an earlier study—. 



n. EXPERIMENTAL DETAILS 

The compound BiPd undergoes polymorphic transfor- 
mation from a-BiPd (monoclinic, P2i)i2 to /3-BiPd 
(orthorhombic,Cmc2i) above 210°C. We have synthesized 
phase pure c-BiPd which has a monoclinic structure with the 
space group P2i with b as its unique axis. Due to its low melt- 
ing point ( 650°C), we have chosen to make the sample us- 
ing a modified Bridgeman technique. The sample was made 
by control heating the individual components (Bi, 99.999% 
pure and Pd, 99.99% pure) in a high purity Alumina crucible 
with a pointed bottom which is kept in a quartz tube that is 
sealed under a vacuum of 1CT 6 mbar. Initially the contents 
were heated upto 650°C (melting point of BiPd) in 12 hrs and 
then kept at 650°C for 12 hrs. Thereafter it was slow cooled 
to 590°C with a rate of l°C/hr and finally the furnace was 
switched off. We obtained high quality poly and single crys- 
tals of few mm size with mass ranging from 10 to 50 mgm. A 
piece of the from the melt was crushed into a fine powder for 
powder x-ray diffraction measurement using Cu Ka radiation 
in a commercial diffractometer. The unit cell of the mono- 
clinic structure (a-BiPd, space group P2i) is shown in Fig. 1. 
The structure consists four inequivalent sites for Bi and four 
inequivalent sites for Pd having 16 atoms in the unit cell. It 
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FIG. 1: (Color online) c-BiPd has Monoclinic(P2|) structure with 16 
atoms in a unit cell (8 formula unit). It contains four inequivalent Bi 
sites and four inequivalent Pd sites. It has no inversion symmetry in 
its monoclinic structure. This structure is stable below 200°C. 



also has alternate layer of Bi and Pd sheets with short Pd-Pd 
distances (shorter than those even in pure Pd metal) with no 
inversion symmetry in the structure. The Rietveld fifcii to the 
powder X-ray data is shown in Fig. 2. The values for the lat- 
tice constants estimated from the fit are a= 5.6284(±0.0004)A, 
b= 10.6615(+0.0004)X, c= 5.6752(±0.0004)A, ar=y=90 and 
y3=101. These values are in agreement with an earlier report- 
except for the value of /?. 

A commercial SQUID magnetometer (MPMS 5, Quantum 
Design, USA) was used to measure the temperature depen- 
dence of the magnetic susceptibility \ m a field of 5 mT for 
temperatures between 1.8 to 20 K to detect the supercondcut- 
ing transition and in a field of 0. 1 T in the temperature range 
from 10 to 300 K. The resistivity was measured using a four- 
probe ac technique on a home built setup with contacts made 
using silver paint on a bar shaped sample 1 mm thick, 10 mm 
long and 2 mm wide. The temperature was measured using a 
calibrated Si diode (Lake Shore Inc., USA) sensor. The sam- 
ple resistance was measured with a LR 700 AC Bridge (Linear 
Research, USA) with a current of 5 mA. The heat capacity was 
measured using a commercial setup (PPMS, Quantum Design, 
USA) in the temperature range from 2 to 200 K. 



FIG. 2: (Color online) Powder X-ray diffraction data of the mono- 
clinic (P2i) ff-BiPd.The solid line is the simulated data using Full- 
Prof(Rietveld program) with the unit cell parameters obtained from 
Ref. 13. 



III. RESULTS AND DISCUSSION 



A. Resistivity studies 

Fig. 3(c) shows the temperature dependence of the resis- 
tivity (p(T)) from 1 .5 to 300 K. The high quality of the sam- 
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FIG. 3: (Color online) Temperature dependence of the electrical re- 
sistivity of the monoclinic (P2i) a-BiPd. The insets a and b are de- 
scribed in the text. The solid lines are fits to the equations described 
in the text. 

pie is clearly evident from the large residual resistivity ratio 
(P300k/P4k) of 160. The inset 3(a) shows the superconduct- 
ing transition below 3.87 K with a width of 10 mK. The inset 
3(b) shows the low temperature resistivity from 4 to 30 K. The 
solid line is a fit to the equation, 



p(T) = po + A XT" 



(1) 



The fit yields a value of 0.3 fi£l cm for po, 9. 2x 10 5 yuf2cm/K 3 
for A and n=3. This T 3 dependence of the resistivity is in- 
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terpreted by the Wilson's theory 15 which takes into account 
the interband s-d phonon induced scattering within the Debye 
approximation. However, the same theory predicts a linear 
temperature dependence for the resistivity at higher tempera- 
tures (T>9d, where Oo is the Debye temperature) which is not 
observed in a-BiPd (see Fig. 3(c)). The Wilson's theory nei- 
ther considers the actual structure of the density of states of 
the electrons at the Fermi level nor the effect of unharmonic- 
ity of the phonon mode. But the contribution from these ef- 
fects in the case of a-BiPd are small. There is yet another 
mechanism (parallel resistor model) suggested by Fisk and 
Webb 16 and later by Wisemann et alii which accounts for the 
significant deviation of the resistivity from the linear tempera- 
ture dependence at high temperatures (100 K<T<300 K). This 
has been seen in many other compounds where the p value is 
rather high (a; 100 pQ. cm). The strong deviation from lin- 
earity and possible tendency towards saturation occur because 
the mean free path becomes short, of the order of few atomic 
spacings. When that happens, the scattering cross section 
will no longer be linear in the scattering perturbation. Since 
the dominant temperature-dependent scattering mechanism is 
electron-phonon interaction here, the p will no longer be pro- 
portional to the mean square atomic displacement, which is 
proportional to T for a harmonic potential. Instead, the resis- 
tance will rise less rapidly than linearly in T and will show 
negative curvature (d 2 p/dT 2 <0). This behavior is also seen in 
our previous studies on silicides and germanidesi^2£. 

Wisemann et describe the p(T) of these compounds 
(which is known as the parallel resistor model) where the ex- 
pression of p(T) is given by, 
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p(T) Pl {T) p II 



(2) 



where p max is the saturation resistivity which is independent of 
temperature and p\ (T) is the ideal temperature-dependent re- 
sistivity. Further, the ideal resistivity is given by the following 
expression, 



Pi CD = po + Ci( 



f )3 f 

Vd Jo 



dx 



(1 - exp(-x))(exp{x) - 1) 



(3) 

where po is the residual resistivity and the second term is due 
to phonon-assisted electron scattering similar to the s-d scat- 
tering in transition metal compounds. 9d is the Debye temper- 
ature and C\ is a numerical constant. Eqn.(2) can be derived 
if we assume that the electron mean free path / is replaced by 
I + a (a being an average interatomic spacing). Such an as- 
sumption is reasonable, since infinitely strong scattering can 
only reduce the electron mean free path to a. Chakraborty 
and Allen 2 ! have made a detailed investigation of the effect 
of strong electron-phonon scattering within the framework of 
the Boltzmann transport equation. They find that the interband 
scattering opens up new nonclassical channels which account 
for the parallel resistor model. The high temperature resistiv- 
ity fit (50 K<T<310 K) shown in the Fig. 3(c) yields a value 
of 164 pQ. cm for p max , 0.33 pQ. cm for po, 76.3 pQ cm/K 3 for 
C p and 168.6 K for do- This value to 9q is close to the value 
obtained from the heat capacity data (described later) which 



suggests that the parallel resistor model can successfully ex- 
plain the high temperature dependence of p(T) of ff-BiPd. 



B. Magnetic susceptibility studies 

The temperature dependence of the magnetic susceptibil- 
ity Cv(T)) of a-BiPd is shown in Fig. 4. The inset shows 
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FIG. 4: (Color online) Temperature dependence of the magnetic sus- 
ceptibility of the monoclinic (P2[) o--BiPd in 1 T. The inset shows 
the temperature dependence of the susceptibility in a field of 5 mT in 
the zero field (ZFC) and field cooled (FC) states. 

the low temperature susceptibility data in 5 mT which re- 
veals the diamagnetic transition of a-BiPd below 3.8 K that 
is in agreement with the resistivity data. The inset also shows 
the zerofield cooled and field cooled susceptibility data which 
elucidates low but significant pinning of the vortices in this 
compound. The high temperature x data of cr-BiPd in its nor- 
mal state is diamagnetic and its value ranges from -2.2xl0~ 5 
emu/mol at 300 K to 1.4xl0~ 5 emu/mol at 1.8 K. The dia- 
magnetic susceptibility at 300 K is indeed surprising given 
that the value of the Sommerfeld coefficient (y=41 mJ/mole 
K 2 ) is quite large. We believe that the observed diamagnetism 
in (T-BiPd is due to the large core diamagnetism of Bi. In gen- 
eral, the observed value of the susceptibility can be written 
as, 



Xobs — Xc< 



XLandau XPauli 



(4) 



Xmre is the core diamagnetism, XLandau is the Landau diamag- 
netism and XPauli is the Pauli paramagnetism. The equation 
(6) can be rewritten as, 



1 m 

Xobs ~ Xcore = XPaid&\ ~ (~ ( )] 

3 m b 



(5) 



Here, XPauli — Na p 2 B S N(Ef ) where is the Avogadro num- 
ber, pb is the Bohr magneton, S is the Stoner factor, N(E/r) 
is the density of states at the Fermi level, m is the free elec- 
tron mass and rm, is the band mass. The estimated value of 
XPauli is 1.3xl0~ 4 emu/mol. If one assumes the core diamag- 
netic susceptibility of Bi as -5xl0~ 5 emu/mol and that of Pd 
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as 0.2x1 -6 emu/mol, we get a value -5.02 xlO~ 5 emu/mol 
as the total contribution to the core diamagnetism in a-BiPd. 
From this we calculate the Stoner enhancement factor S as 3. 
A weak temperature dependence of x(T) could arise due to 
two reasons. One of them could be the presence of magnetic 
impurities at the 100 ppm level in the compound. This could 
account for the observed small temperature dependence of the 
X(T). However, we believe that impurity concentration of this 
amount either in Bi or in Pd is unlikely (given that the limit 
of the impurities both in Bi and Pd elements imposed by the 
manufacturer is <10 ppm). The second reason could be the 
temperature variation of the density of states at the Fermi level 
which results in a temperature-dependent Pauli spin suscepti- 
bility as seen in some of the A-15 compounds. Knight-shift 
measurements would be useful to resolve this issue. 



C. Heat capacity studies 




FIG. 5: (Color online) Temperature dependence of the heat capacity 
of the monoclinic (P2i) o--BiPd from 2 to 35 K. The insets are ex- 
plained in the text. Solid line is a fit to the equation(6) described in 
the text. 

The temperature dependence of the heat-capacity (C p ) from 
2 to 30 K of or-BiPd is shown in Fig. 5. The inset shows the 
low temperature C p vs T data. The jump in C p at 3.7 K (AC= 
30 mJ/mol K) clearly shows bulk superconducting ordering in 
this sample below this temperature. The temperature depen- 
dence of C p is fitted to the expression, 



C p = yT + j3 r j 



(6) 



where y is due to the electronic contribution and /3 is due to 
the lattice contribution. The value of the ratio AC p /yT c is 0.2 
which is significantly reduced from the BCS value of 1.43. 
Low values of AC p /yT c have been observed before in the heat- 
capacity study of several superconducting compound o 22,23 . 
According to these studies the reduced jump across T c could 
arise from extrinsic effect (such as inhomogeneity in the sam- 
ple or magnetic impurities) or from intrinsic effect (such as 



the existence of regions which do not participate in supercon- 
ductivity). In our sample of ff-BiPd, we estimate the impurity 
content to be less than 0.1 % by volume and the sharpness 
of the superconducting transition also suggests good homo- 
geneity.lt is possible that the reduced jump could arise from 
two-band superconductivity where one band remains normal. 
However, detailed Fermi surface measurements are required 
before we can analyze the data in terms of this model. The fit 
to the heat capacity data using the eqn. 6 in the temperature 
range from 5 to 20 K yielded 41 mJ/mol K 2 and 0.9 mJ/mol 
K 4 for y and /?, respectively. The y value has been obtained 
by matching the entropy of normal and superconducting state 
at T c as suggested by Stewart et a£&. From the /3 value of 0.9 
mJ/mol K 4 , we estimate the 6q to be 162 K using the relation, 



12 n 4 N A n k B 
5^ 



1/3 



(7) 



where is the Avogadro's number, n is the number of atoms 
per formula unit, and kg is the Boltzmann's constant. 



1. Upper critical field studies 

The estimation of the upper critical field (H e 2) value at a given 
temperature has been made by measuring the resistance of the 
sample under a given magnetic field. The transition tempera- 
ture in a given field is defined as the temperature which corre- 
sponds to the mid point of the resistance jump. The tempera- 
ture dependence of H C 2 is shown in Fig. 6. It is well known 
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FIG. 6: (Color online) Temperature dependence of the upper criti- 
cal field (H„(T)) of a-BiPd. The estimation was made using both 
resistivity and magnetization measurements as described in the text. 

that in non-magnetic superconductors, the magnetic field in- 
teracts with the conduction electrons basically through two 
different mechanisms. Both lead to pair breaking and even- 
tually destroys the superconductivity at a given field which 
is known as the critical field. One of this mechanisms arises 
due to the interaction of the field with the orbital motion of 
the electrons (orbital pair breaking) and the other is due to the 
interaction of the field with the electronic spin (Pauli param- 
agnetic limiting effects). Orbital pair breaking is the dominant 
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mechanism at low fields and at very high fields, Pauli param- 
agnetic effect limits the upper critical field. We have fitted this 
temperature dependence of H C 2 to the theory of clean type-II 
superconductor^ 2 - 7 -. We obtain a value of 0.16 for A, r , 0.69 T 
for H c ,(0) and 3 kOe/K for dH f 2/dT near T c . The value of the 
Pauli paramagnetic limiting field (H/> a »/;=L84xT c in (Tesla)) 
for a-BiPd is very large (7.05 T) compared to the estimated 
value of the upper critical field at K. From the GLAG 2 ^— 
theory, we get the value of the upper critical field at K as, 



H C2 (0) = 0.72 (dHJdT)T e 



(10) 



From this we get a value of 0.8 T for H Cl (0) which is larger 
than the extrapolated experimental value of H c , . However, 
assuming the clean limit for the type-II superconductor, one 
can also estimate the dH t 2/dT using the relation, 



dH c JdT 



9.55 10 24 y 2 



T c [« 2/3 (S/S F )T 2 (in Oe/K) (11) 



where y is the electronic heat capacity coefficient 
(ergs/cm 3 K 2 ), n is the conduction electron density in 
units of crrT 3 and S/Sy is the ratio of the area of the Fermi 
surface to that of a free electron gas of density n. Substituting 
the values of y,n and assuming a simple model of spherical 
Fermi surface (S/S F =1), we get a value 2.08 kOe/K which 
is larger the value (3 kOe/K) obtained from the experiment. 
The reason for this large discrepancy is not understood at this 
moment though similar anomalies in the value of dH C 2/dT 
have been reported in earlier studie o 31,32 . In those earlier 
reports they have studied strong-coupled superconductors 
which have complex phonon spectra. In that case, utilizing 
the theory to analyze is not strictly valid as the theory assumes 
electron interact via weak-coupling BCS type interaction 
potential and have spherical Fermi surface. 



D. Estimation of normal and superconducting states 
parameters 

Using the value of 9d and T c , we can estimate the electron- 
phonon scattering parameter, A, from the McMillan's theory 33 
where A is given by, 



1.04 + fi* MfipjlAS T c ) 
(1 - 0.62 fi*) ln(8 D /l.45 T c ) - 1.04 



(12) 



Assuming //=0.1, we find the value of A to be 0.72 which 
puts a-BiPd as an intermediate coupling superconductor. On 
the basis of purely thermodynamical arguments, the thermo- 
dynamic critical field at T=0 K (Hc(0)) can be determined by 
integrating the entropy difference between the superconduct- 
ing and normal states. From our experimental heat capacity 
data, we obtain a value of 2000 Oe for Hc(0). One can also 
calculate the thermodynamical critical field Hc(0) from the 
expression 2 ^, 



H c (0) = 4.23 y l/2 T, 



(13) 



where y is the heat capacity coefficient(erg/cm 3 K 2 ). This 
gives a value of Hc(0) as 2050 Oe. We can estimate the 



Ginzburg-Landau coherence length at T=0 K from the re- 
lation, 



f<a(0) = 



5.87 1Q- 17 » 2/3 (575V) 
J To 



(14) 



where y is electronic heat capacity coefficient (erg/ cm 3 K 2 ) 
and S/Sf is the ratio of the area of the Fermi surface to that 
of a free electron gas of density n. This equation yields a 
value of 205 X for ^ci(O). If one uses the standard rela- 
tion, H c ^(0) — t „ fabo? where (f>o is the flux quantum 

(2.07xl0~ 7 oersted/cm 2 , we get a value of 218 Xfor £ GL (0) 
which is closer to the estimated value from the equation (14). 
Using this value for ^cl(O), one can estimate k(0) as 2.44 
(since k(0) = 5" (Q) ). However, from the GLAG theory for 

clean superconductors 2 ^, we know that k(Q) = "^^^jgW^ 
(where k(0) - A cl (0)/^ C l(0). Using that, we get the k(0) value 
as 2.7 which is close to the value estimated earlier. From 
the value of ^cl(O) = 218 A, we get a value of 532 A for 
the Ginzburg-Landau penetration depth at K (/Igl(0)). The 
lower critical value can be determined by using the relation, 



HciQS) 



Hc (0) ln[ K (0)] 
2 1 / 2 k(0) 



(15) 



which yields a value of 516 Oe for the lower critical field at 
K. This value of Hci(0)is in agreement with magnetization 
measurements on the same sample. The enhanced density of 
states can be calculated using the relation, 



N*(E F ) = 0.2121 y/N 



(16) 



where N is the number of atoms per formula unit and 
y is expressed in mJ/mol K 2 . The value of N*(E F ) 
is 4.2 states/(eV atom spin-direction) and the value of 
the bare density of states N(E F ) = N*(E F )/(l+/l) = 2.4 
states/(eV atom spin-direction). 

The parameters are calculated using the Ginzburg-Landau 
theory for clean type-II superconductors. To verify the self 
consistency of our approach, we have estimated the mean free 
path (/) of our sample using the expression, 



I = 1.27 10 4 [p« 2/3 (S/S F )Y 



(17), 



where n is the conduction electron density in units of crrT 3 
and S/Sf is the ratio of the area of the Fermi surface to that of 
a free electron gas of density n. If one assumes a simple model 
of spherical Fermi surface (S/Sj?=l), the value of / would be 
1965 X and one can also calculate the value of the BCS co- 
herence length (for S/Sf=l) from the expression, 



ft 



7.95 10~ 17 [n 



2/3 



}{yT c y 



(18) 



where y is expressed in ergs/cm 3 K 2 . The value of £o was 
found to be 277 A which is much lower than I which implies 
ff-BiPd is a clean type II superconductor. Moreover, we get 
the value of 0.13 for A tr using the equation A, r = 0,882 — , 
which is close to the value of 0.16 obtained from the HW fit 
earlier. 
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E. Conclusion 

Noncentrosymmetry (NCS) introduces an electrical field 
gradient in the crystal which creates a Rashba-type antisym- 
metric spin-orbit coupling 2 . An inherent feature is then the 
mixing of spin-singlet and spin-triplet Cooper pairing chan- 
nels which are otherwise distinguished by parity. This mixing 
of pairing states is expected to cause a two-component order 
parameter. Such a behavior is observed in heavy fermion com- 
pounds like CePt 3 Sii, CeRhSiA CelrSi^ and Ulri. Apart 
from heavy fermion properties, the occurrence of supercon- 
ductivity in materials without inversion symmetry deserves 
its own merit. Therefore, we aim to find NCS superconduct- 
ing materials without strong correlations among d or f elec- 
trons. This may set the stage for reliable electronic struc- 
ture calculations, proving the splitting of bands due to miss- 
ing inversion symmetry. Previous examples of superconduc- 
tors without strong electron correlation, where, however, NCS 
was not a central issue of investigation, are binary R2C3 with 
R=La or or Cd 2 Re 2 7 with Tc=l K 7 . More recently 



the solid solution Li2Pd,Pt3B was demonstrated to show 
superconductivity . While l^Pd^B is accounted for in terms 
of a conventional BCS-like superconductor, Li2Pt3B refers to 
unconventionally^ presumably due to a very large spin-orbit 
coupling. A recent study on Li2Pd.vPt3_.vB by muon-spin rota- 
tion and specific heat, however, suggests that the whole series 
comprises single-gap s-wave superconductivity 34 . As far as 
ff-BiPd is concerned, all experimental parameters and anal- 
yses point to a BCS-like superconductivity, which is based 
on spin-singlet pairing with a fully gapped DOS at the Fermi 
level. Therefore, a-BiPd is classified as a type II supercon- 
ductor in the clean limit. Here, the absence of inversion in the 
crystal structure does not give rise to anomalous features of 
the superconducting state such as a mixing of spin-singlet and 
spin-triplet pairs, as well as nodes in the superconducting gap 
at the Fermi surface. The overall effect of non-inversion sym- 
metry seems to be of minor importance for the superconduct- 
ing properties. In this respect, c-BiPd seems to be in line with 
conclusions drawn for Li2Pd3B and BaPtSi^. Band structure 
calculations of c-BiPd are in progress. 
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